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Abstract 

We study a possibly integrable model of abelian gauge fields on a two-dimensional 
surface M , with volume form ^ . It has the same phase space as ideal hydrodynam- 
ics, a coadjoint orbit of the volume-preserving diffeomorphism group of M . Gauge 
field Poisson brackets differ from the Heisenberg algebra, but are reminiscent of 
Yang-Mills theory on a null surface. Enstrophy invariants are Casimirs of the Pois- 
son algebra of gauge invariant observables. Some symplectic leaves of the Poisson 
manifold are identified. The Hamiltonian is a magnetic energy, similar to that of 
electrodynamics, and depends on a metric whose volume element is not a multiple of 
^. The magnetic field evolves by a quadratically non-linear 'Euler' equation, which 
may also be regarded as describing geodesic flow on SDiff (M, /i) . Static solutions 
are obtained. For uniform /i, an infinite sequence of local conserved charges begin- 
ning with the hamiltonian are found. The charges are shown to be in involution, 
suggesting integrability. Besides being a theory of a novel kind of ideal flow, this is 
a toy-model for Yang-Mills theory and matrix field theories, whose gauge-invariant 
phase space is conjectured to be a coadjoint orbit of the diffeomorphism group of a 
non-commutative space. 
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1 Introduction and Summary 



The classical theory of gauge fields"'^ Ai{x^ ,x'^,t) and A2{x^ ,x'^ ,t) we study in this paper, may 
be summarized in four equations. The hamiltonian is a gauge-invariant magnetic energy, 
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where B = diA2 — i92^i is the magnetic field, p is a given volume element and gij a fixed 
metric on a two dimensional surface M such that a = p^/g is not a constant {g = det gij). The 
Poisson bracket between gauge fields is a Lie algebra (independent of the metric) 



{A{x),A,{y)} = S\x-y) 



-1 
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Hamilton's equation for time evolution of gauge-invariant observables, / 
the magnetic field evolves according to a non-linear 'Euler' equation 

B = V{B/p) X V{B(t/p). 



{H, /} implies that 
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Viewed as a rigid body for the group of volume preserving diffeomorphisms of M , the inverse 
of the inertia tensor is a twisted version of the Laplace operator 
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This theory has the same phase space and Poisson brackets as 2 + 1 ideal(inviscid and volume 
preserving) hydrodynamics, but a different hamiltonian. It may be integrable, since we find an 
infinite number of conserved quantities Hn = J{B/p)"' a p cPx as well as an infinite number 
of Casimirs In = J{B/p)^ p cPx for uniform p. It is remarkable that one can make this 
modification to 2 + 1 ideal fiow, which is sometimes studied as a toy-model for turbulence, to 
get a potentially integrable system. However, our original motivation for studying this model 
was different. We argue below that it is the simplest 'gauge theory' that shares some quite deep, 
though unfamiliar, mathematical features of Yang- Mills theory. 

The formulation of Yang-Mills theory in terms of gauge-invariant observables, and the de- 
velopment of methods for its solution are important and challenging problems of theoretical 
physics, since all the experimentally observed asymptotic states of the strong interactions are 
color-singlets. This problem has a long history stretching at least as far back as the work of 
Mandelstamj^. Wilson loops are a natural choice for gauge-invariant variables, but they have 
trivial Poisson brackets on a spatial initial value surface, since the gauge field is canonically 
conjugate to the electric field on such a surface. More recently, it has been shown by Rajeev 
and Turgut [21 El , that Wilson Loops of 3 -|- 1 dimensional Yang- Mills theory on a null initial 
value hypersurface satisfy a quadratic Poisson algebra with no need for electric field insertions. 
This is because the transverse components of the gauge field satisfy a non-trivial Poisson alge- 
bra among themselves, as opposed to the situation on a spatial surface. The Poisson algebra 
of Wilson loops is degenerate due to Mandelstam-like constraints. The gauge invariant phase 
space of Yang-Mills theory is conjectured to be a coadjoint orbit of this Poisson algebra. It is 

^After gauge-fixing, there will be a single propagating field degree of freedom. 
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still a challenge to write the hamiltonian in terms of these variables. However, this has been 
possible in dimensionally reduced versions^ such as adjoint scalar field theories coupled to quarks 
in 1 + 1 dimensions, as shown by Lee and RajeevO IHI- In conjunction with 't Hooft's large 
approximation 0, viewed as an alternative classical limit, this is an approach to better under- 
stand the non-perturbative dynamics, especially of non-supersymmetric gauge theories. In such 
an approach to 1 + 1 QCD, the phase space of gauge-invariant meson variables is an infinite 
Grassmannian, a coadjoint orbit of an infinite dimensional unitary group [Hj. This allows one to 
understand baryons as well as mesons in the large limit, going beyond the early work of 't 

Hooftpiaiiniiii]. 

However, the groups and Lie algebras whose coadjoint orbits are relevant to matrix field the- 
ories and Yang-Mills theory are poorly understood non-commutative versions of diffeomorphism 
groups'^. In the case of a multi-matrix model, the group is, roughly speaking, an automorphism 
group of a tensor algebra. The Lie algebra is a Cuntz-type algebra which can be thought of 
as an algebra of vector fields on a non-commutative space ^1 El El El However, it is still 
very challenging to find the proper mathematical framework for these theories and develop ap- 
proximation methods to solve them even in the large limit. To develop the necessary tools, 
it becomes worth while to practice on simpler theories whose gauge invariant phase space is 
the coadjoint orbit of a less formidable group. Here, we take a step in this direction by study- 
ing an abelian gauge theory whose phase space is a coadjoint orbit of the volume preserving 
diffeomorphism group of a two dimensional surface. 

To put these remarks in perspective, recall the common classical formulation of Eulerian 
rigid body dynamics, ideal hydrodynamics, the KdV eauation P^ [TBI [T7j and the large limit 
of two dimensional QCDjS]. The phase space of each of these theories is a symplectic leaf of a 
degenerate Poisson manifold, which is the dual G* of a Lie algebra. G* always carries a natural 
Poisson structure. Symplectic leaves are coadjoint orbits of a group G acting on the dual of 
its Lie algebra G* . On any leaf, the symplectic structure is given by the Kirillov form. The 
appropriate groups in these examples are 50(3) , the volume preserving diffeomorphism group of 
the manifold upon which the fluid flows, and the central extensions of Diff (5"^) and of an infinite 
dimensional unitary group, respectively. The coadjoint orbits for the rigid body and 2d QCD are 
well-known symplectic manifolds: concentric spheres and the infinite dimensional Grassmannian 
manifold. The observables in each case are real-valued functions on G* ■ The Poisson algebra 
of observables is degenerate, i.e. has a center consisting of Casimirs. The symplectic leaves 
can also be characterized as the level sets of a complete set of Casimirs. In each case, the 
hamiltonian is a quadratic function on the phase space and classical time evolution is given 
by Hamilton's equations. Hamilton's equations are non-linear despite a quadratic hamiltonian, 
since the Poisson brackets of observables are more complicated than the Heisenberg algebra. In 
exceptional cases such as the rigid body and the KdV equation, these non-linear equations are 
exactly integrable. In other cases, it is useful to develop approximation methods to solve them, 
that are adapted to the geometry of the phase space. 

Our earlier remarks indicate that it may be fruitful to regard Yang-Mills theory and matrix 
field theories as Hamiltonian dynamical systems along the lines of the more well known ones 
listed in the last paragraph. As a toy-model in this direction, we seek a gauge theory where 

^See also the work of Karabali, Nair and Kim who have made significant progress with a gauge invariant 
Hamiltonian approach to 2 + 1 Yang-Mills theory 

^This is not the structure group (sometimes called the gauge group) of the theory, which is still SU{N) or 
U {N) . The gauge group plays little role in a gauge-invariant formulation of the theory. 
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the gauge fields satisfy a closed Poisson algebra, without any need for electric fields. We want 
a theory whose phase space is a coadjoint orbit of an ordinary difFeomorphism group, which 
is simpler than its non-commutative cousins. We would also like to understand in more detail 
the structure of the Poisson algebra of gauge-invariant observables, work out the equations of 
motion and try to solve them. 

In this paper, we identify a classical theory of abelian gauge fields in two spatial dimensions, 
different from Maxwell theory. In particular, it is not Lorentz covariant, indeed, time plays 
the same role as in Newtonian relativity. The theory is defined by a two dimensional manifold 
M, a volume form /i and a metric gij whose volume element is not a multiple of //. The 
hamiltonian is a gauge- invariant magnetic energy, much like that of Maxwell theory. Unlike in 
electrodynamics, the gauge field is a 1-form on space, rather than on space-time. Thus, even 
before any gauge fixing, the gauge field has no time component. There is a magnetic field B , but 
no electric field, so to speak. After gauge fixing, there remains only one dynamical component 
of the gauge field. In this sense, the theory has the same number of degrees of freedom as 2 -|- 1 
electrodynamics. However, though the hamiltonian is quadratic in the gauge fields, the classical 
theory is nonlinear due to the 'non-canonical' Poisson algebra of gauge fields. Equations of 
motion are non-linear and comparable to those of a 2-1-1 dimensional non-abelian gauge theory 
or ideal hydrodynamics. 

The phase space of the theory is a coadjoint orbit of the volume preserving diffeomorphism 
group SDiff(M, /x) of the spatial two-dimensional manifold. Roughly speaking, this means that 
SDiff (M, /i) is a symmetry group of the Poisson algebra of observables. The gauge group (struc- 
ture group) of the theory is U{1). The inspiration for this lies in ideal hydrodynamics. Indeed, 
even before the diffeomorphism group of a manifold appeared in general relativity, it was relevant 
as the configuration space of a fiuid. The theory we study is not the same as, but is motivated by 
2-|- 1 dimensional ideal hydrodynamics, regarded as a hamiltonian svstem [TFl 1191 ITCl I2U1 OTl I22j . 
Though we arrived at it as a toy model for Yang-Mills theory, it turns out to have a nice geomet- 
ric and possibly even integrable structure. We find two infinite sequences of conserved charges. 
The first set are Casimirs, analogues of the enstrophy invariants of ideal hydrodynamics. In 
addition, we find another infinite set of conserved charges which are not Casimirs but are in 
involution. The theory we study here can also be regarded as a theory of geodesies of a right- 
invariant metric on the volume preserving diffeomorphism group of a two dimensional manifold. 
However, the right invariant metric on SDiff (M,/x) implied by our hamiltonian is different from 
that arising in ideal hydrodynamics (the metric leading to ideal Euler fiow) as well as the 

metric leading to averaged Euler fiow 23]. 

Another way to view the current work is to recall that adding supersymmetry usually gives 
greater analytical control over gauge theories. But there may be other modifications of gauge 
theories that also lead to interesting toy-models or enhanced solvability. Our investigation 
concerns one such novel modification of gauge field Poisson brackets. 

In Sec. [21 we introduce the space of abelian gauge fields Aidx^ on a two dimensional surface 
M as the dual of the Lie algebra SVect(M, ^) of vector fields preserving a volume element 
fj- = p cPx. This 'duality' is known in hydro dynamics ^3]. The differentials d/* = P~^jj- of 
differentiable gauge-invariant observables f{A) are shown to be volume preserving vector fields. 
In Sec. 01 we give the Poisson structure on gauge-invariant observables 

{/..} . /.v^.[.-^s.(.-'^)-p-^a,(.-'i|-)] (5) 
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turning the space of gauge fields into a Poisson manifold. The Poisson brackets of gauge fields 
are obtained explicitly H35() and compared with those of Yang-Mills theory on a spatial and null 
initial value hyper surf ace. 

In Secl^we give the coadjoint action of SDiff(M, /x) and its Lie algebra SVect(M, ^u) on the 
Poisson manifold of gauge fields SVect(M, /x)* , and show that the action is canonical i.e. pre- 
serves the Poisson structure. The moment maps generate the coadjoint action. The symplectic 
leaves of the Poisson manifold are coadjoint orbits. The enstrophy invariants of hydrodynamics 
In = J]^.j{dA/ fi)'^ fj, are an infinite sequence of Casimirs of the Poisson algebra. The coadjoint 
orbits of closed gauge field 1 -forms are shown to be finite dimensional. Single-point orbits for 
simply connected M are found. We argue that all other orbits are infinite dimensional and try 
to characterize their isotropy subalgebras as well as tangent spaces. 

In Sec. El we first review the choice of hamiltonian leading to ideal Eulerian hydrodynamics 
in 2+1 dimensions. Then we propose a different gauge-invariant hamiltonian depending on both 
fj, and a metric gij , by analogy with the magnetic energy of Maxwell theory. 

where F = dA is the field strength, *F is its Hodge dual, and J7g is the volume element of the 
metric gij . Here a = (///Og)^ = p^/g, g = det gij and p is the density associated to fi. H is 
shown to determine a non-negative inner product on the dual of the Lie algebra SVect(-/Vf, //)* 
and an inverse 'inertia tensor' by analogy with the rigid body. If M is simply connected, the 
inverse inertia operator is non-degenerate and could be inverted to get an inner product on the 
Lie algebra SVect(M, /i). This could be extended to the diffeomorphism group SDiff(M, /x) by 
right translations. Thus, the magnetic energy should define geodesic flow on SDiff(M, /i) with 
respect to a right-invariant metric different from that coming from Eulerian hydrodynamics. 

In Sec. iniwe find the equation of motion for the magnetic field B = e^^diAj , B = 'V{B/p) x 
V{Bp/g). This simple quadratically non-linear evolution equation is strikingly similar to the 
Euler equation of a rigid body L = L x Q, L = lil.. It can be regarded as the 'Euler equation' 
for the group SDiff(M, /x) with hamiltonian given above. Remarkably, for a uniform measure 
p, we find an infinite sequence Hn of conserved charges in involution, which are not Casimirs. 
The hamiltonian is ^ . 

Hn= [ {dA/p^ap, n = 1,2,3,... (7) 

In Sec. [7|we find some static solutions of the equations of motion. We show that for circularly 
symmetric p and g , every circularly symmetric magnetic field is a static solution. We generalize 
this to the non-symmetric case as well. We also find a one parameter family of static solutions 
that are local extrema of energy even with respect to variations that are not restricted to the 
symplectic leaf on which the extremum lies. Some ideas for further study are given in Sec. |HJ 

Volume preserving diffeomorphisms and gauge theories have appeared together previously 
in the literature (see for example Ref. |241 125j ). Our investigation seems quite different, since 
SDiff (M, p) is not the gauge group of our theory but rather a symmetry of the Poisson algebra. 
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2 Volume Preserving Vector Fields to Gauge Invariant Observ- 
ables 



2.1 Lie Algebra of Volume Preserving Vector Fields 

Let M be a surface with local coordinates , to be thought of as the space on which a fluid 
flows. A vector field on M is regarded as the velocity field of a fluid at a particular time. The 
space of all vector fields on M forms a Lie algebra Vect(M) with Lie bracket 

[u, vY = u^djv' - v^dju' (8) 

Vect(M) is the Lie algebra of the diffeomorphism group Diff(M). Conservation of the mass of 
the fluid during its flow implies the continuity equation for its density p{x, t) 

^ + V.(pu) = (9) 

We are interested in flows where the density at any point of space does not depend on time. 
Using the continuity equation, this becomes V- (pu) = 0. We call such a flow volume preserving. 
Geometrically, we are considering a flow that generates diffeomorphisms of M that preserve a 
given volume form^ jC^p = 0. The density is constant along integral curves of u. To see the 
equivalence of this with the continuity equation for a volume preserving flow, recall that 

CuH = {diu + iud)p = d{iulj) (10) 

where z„ is the contraction with u. Here dp = since p is a volume form. In local coordinates 
p = ^pijdx^ A dx^ where pjj = — Pji = eijp and is antisymmetric with €12 = 1. So 
p = p(x) d'^x where dx^ A dx"^ = (fx. Then i„p = ^pij{u''dx^ — dx'^u^) = piju^dx^ , so that 

CuP = d{iup) = dk{iiijU^)dx^ A dx-^ = di{pu^)dx^ A dx^ (11) 

Thus >Cup = becomes V • (pu) = 0. We will use the terms volume preserving and area 
preserving interchangeably since M is a two dimensional surface. Some of what we say has a 
generalization to higher (especially even) dimensional M . 

The properties Cau+f3v = o:Cu + and = CuCv — C-uCu ensure that the space of 

volume-preserving vector fields Q = SVect (M,p) forms a Lie subalgebra of Vect(M). It is the 
Lie algebra of the group of volume preserving diffeomorphisms G = SDiff(M, p). 

Volume preserving flow is a mathematical implementation of the physical concept of incom- 
pressible flow, which occurs when the fluid speed is small compared to the speed of sound^. In 
particular, shock waves cannot form in incompressible flow, since shock waves involve super- 
sonic flow. Under ordinary conditions, air flow in the atmosphere is incompressible. Volume 
preserving flow is sometimes referred to as divergence-free flow. 

*A volume form must be non-degenerate. In 2 dimensions it is the same as an area form or a symplectic form. 

^Some authors consider only the special case where density is a constant, V • m = 0. Note also that the same 
fluid may support both compressible and incompressible flow under different conditions, so our definitions refer 
to the flow and not just to the fluid. 
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If M is simply connected, the volume preserving condition di{pu^) = may be solved in 
terms of a stream function ip satisfying = dip . \s a, scalar function on M that serves as 
a 'potential' for the velocity field. In local coordinates 



iuP- = dip ^ PijU^ = djip (12) 
Since // is non-degenerate (p 7^ 0), it can be inverted p~^£^^ p^jk = —^k where e'-' is a constant 



antisymmetric tensor with e^-^ = 1, e^^Cjk = —^l- This does not require a metric on M. Then 



u' = p-^e'^djip. (13) 

u determines ip up to an additive constant, which can be fixed by a boundary condition. If M 
is simply connected, then SVect(M, /i) may be identified with the space of stream functions. 
Suppose two volume preserving vector fields u, v have stream functions ipu and ip^ , 

u' = p-h'^djiPu, v' = p-h'^djiPu (14) 

Then their Lie bracket \u, v\ has stream function p^^Vipv x Vipu '■ 

[u,vY = p-'e''di(^p~'e^\djiP,){dM 

i^iuM = p-'e^HdjA){dM ' (15) 



2.2 Abelian Gauge Fields as the Dual of SVect(M,/i) 

The Lie algebra SVect(M, /i) is akin to the Lie algebra of angular velocities of a rigid body. The 
angular momenta are in the dual space to angular velocities, and satisfy the angular momentum 
Poisson algebra. As explained in0 the dual of any Lie algebra is a Poisson manifold. This is 
interesting because the observables of a classical dynamical system are real-valued functions on 
a Poisson manifold. The dual of the Lie algebra Q = SVect(M, /i) is the space of abelian gauge 
fields modulo gauge transformations, 

g* = SVect(M,^)* = n^{M)/dn°{M). (16) 

This fact is well-known in hydrodynamics (see Ref. JH]), though it is usually not thought of in 
terms of gauge fields. To see this duality, we define the pairing (A, u) between gauge fields and 
volume preserving vector fields by integrating the scalar A{u) with respect to p 

{A,u) = puiA) = [ A{u)p = [ Aiu' p d^x (17) 

The pairing pu{A) is also called the moment map. It is a gauge-invariant pairing. Under a 
gauge transformation A^ A' = A + dA for any scalar A(x) 

Pu{A') - pu{A) = f (a,A)nV = - / Adi{pu')d^x = (18) 
Jm J 

since u is volume preserving. We assume that gauge fields and gauge transformations A vanish 
on the boundary dM or at infinity. We make no such assumption about the vector fields. 

Gauge Fixing: It is occasionally convenient to 'gauge-fix', i.e. pick a coset representative 
for (M) / dCl^ (M) . Under a gauge transformation, A[ = Ai + diA. We can pick A such that 
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A[ = Ai + diA = , so that we are left with only one component of the gauge field A2 . We can 
still make an -independent 'residual' gauge transformation, A2 = A'2 + d2A[x'^) to eliminate 
any additive term in A2 depending on x'^ alone. Suppose we have gauge fixed on a particular 
spatial initial value surface at time t = 0. Unlike in Yang-Mills theory, the equations of motion 
of our theory are purely dynamical. They only evolve the gauge-fixed fields forward in time, 
and do not contain any further constraints. In effect, after gauge fixing, we will be left with one 
propagating field degree of freedom. 



2.3 Differentials of Gauge Invariant Charges are Volume Preserving Vector 
Fields 

We should regard SVect(M, /x)* , the space of gauge fields modulo gauge transformations, as the 
Poisson manifold of some dynamical system. Real-valued functions on this space (i.e. gauge- 
invariant functions f{A)) are the observables. Given such an f{A), we can define its differential 
df{A) 

(df{A)y=p-Hx)j^^^p-'d^f (19) 

For each equivalence class of gauge fields [A] = {A\A ~ A+dA} , the differential^ defines a vector 
field df^di on M . If f{A) is non-linear, the vector field {df{A)y(x) changes as A £ Q* changes. 
Suppose f{A) is gauge invariant and differentiable. Then we can show that its differential d/* 
is a volume preserving vector field on M: di{p dP) = 0. To see this, note that gauge invariance 
implies that the change in / under any gauge transformation 5Ai = diA must vanish 

Since A(x) is arbitrary, it must follow that di{5f/5Ai) = 0. So the differential of a gauge 
invariant function can be regarded as an element of the Lie algebra Q = SVect(M, p) . 

The simplest gauge- invariant observable is the field strength 2 -form 

F = dA = ^Fijdx' A dx^; Fij = diAj - djAi = eijB (21) 
where B = e^WiAj = diA2 — 82 Ai is the magnetic field. The differential of F 

{{dF){A))>' = P-'^^^ = P-' {spAx -y)- 5^d,5{x - y)) (22) 
is a volume preserving vector field on M for each A: 

dk{p{dFf) = dk (^<5}-a,5(x -y)- 6fd,6{x - y)^ = {did, - d,di)6{x -y) = 0. (23) 

®The differential {df{A)y is regarded as a vector field on M for each A and should not be confused with the 
closely related exterior derivative df , which is a 1-form on Q* . However, as we will see later (SecHig, on any 
symplectic leaf of SVect(M, /i)* with symplectic form u; , the 1-form df determines the canonical vector field Vf 
via u>{Vf, .) = df{.) . V/ is a vector field on the leaf, and its relation to the differential is Vf{A) — ad'^^A. 
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Similarly, the differential of the magnetic field 

iidBXA))'^ = p-'^j^, = p-'e''d,S^S{x -y)= p-'e'^^d^Six - y) (24) 

is volume preserving dk{p{dB)^) = e^^ did]^5{x — y) = 0. We can regard the moment maps /Xu(^) 
as linear gauge- invariant observables. The differential of Hui-A) is the volume preserving vector 
field n, for all gauge fields A. 

Other gauge invariant observables f{A) we will be interested in are 'charges': integrals 
over M with respect to /x, of a local gauge- invariant scalar function T . T can depend on A 
only through the field strength two form F = dA. The analogue of the Chern-Simons 3-form, 

vanishes identically since A is a 1-form on space, not space-time. The quotient of dA and the 
non-degenerate volume 2-form {dA/ p) is a scalar function on M . Then 

f{A) = J T{a, {dA/p), v%idA/p), w'^didj{dA/p), . . .) p (25) 

where cr is a scalar function and v^,w^^ etc arc arbitrary but fixed contravariant tensor fields. 
We can get an explicit formula for the differential of such a gauge invariant charge. Using 
dA/p = B/p and 

^^^._..a.(y(.-,))... (.0) 

we get upon integrating by parts, 



df{Ay = p-'^ = p-'s'^dj 



SA. 



{d{B/p)) {p^^iddkiB/p)))^--- 



(27) 



Due to the anti-symmetry of e*-' , it follows that df is volume preserving di(p dP) = 0. Moreover, 
if / is gauge-invariant and of the form assumed above, then its differential is also gauge-invariant. 

Two families of gauge invariant charges which play an important role in our theory are 
and Hn defined below. Let 

In{A) = f {dA/pr p= f{B/pr pd^x= f f-Y \a, n = 1,2,3,... 

Jm J JM \ p J 

^"^^-y = = ^^^''^^■((^/'^)"~') (28) 

P O^i p 

Their differentials are volume preserving di{p{dIn{A)y) = since e*-' is antisymmetric. Note 
that Ii = Jj^ dA = . We assume B vanishes sufficiently fast at infinity and do not consider 
/„ for n < 0. Given a scalar function cr on M we can construct additional gauge- invariant 
charges. These are similar to the /„ , except that we multiply by a before integrating over M 

HniA) = I {dA/pTap = I {B/pT a p d\ (29) 

More generally, one can replace {B / pf' by an arbitrary function of Bj p. The differential of 
is volume preserving: 

{dHr,{A)y = — dMB/pr~')- (30) 
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Gauge-invariant observables depending on the volume form and a metric gij on M will play 
an important role in determining the dynamics of our theory. They are given in Sec. 15.21 

There are other interesting gauge invariant observables such as the circulation C^{A) = 

Jq Ai '^'^Jf^ ds and its exponential, the Wilson loop. Since these observables are concentrated 
on one dimensional curves on M, they may fail to be differentiable and their differentials exist 
only as distributional vector fields. They require a more careful analysis and are not considered 
in this paper. Henceforth, when we say observable, we will mean gauge invariant observables 
that are differentiable, in which case, their differentials are guaranteed to be volume preserving 
vector fields on M . 



3 Poisson Brackets 

3.1 Definition of Poisson Bracket on g* = SVect(M,/i)* 

The Lie algebra structure of volume preserving vector fields Q = SVect(M, /x) can be used 
to define a Poisson structure on its dual (see Ref. [211 l2Zl 1^)- The dual space G* = 
SYect{M, n)* = il,^ (AI) / dO,^ (M) of gauge field 1-forms then becomes a Poisson manifold. Ob- 
servables are real-valued functions on it. The Poisson bracket (p.b.) between gauge-invariant 
observables f{A) and g{A) with volume preserving differentials df and dg, is defined using the 
pairing (^4, u) between Q and G* and the Lie algebra bracket [df, dg] : 

{/, g}{A) ^ {A, [df, dg]) = j Ai[df, dg]' ^jl = j d^xpAi[df^ djdg' - dg^d.df'] 

- J dxpA^p ^^^d,^p ^^J p ^ 

{Pf)d,{p~'6'g)-{Pg)d,{p-'6'f) 



Ai 



d^x (31) 



where (5* = . The antisymmetry, linearity and Jacobi identity follow from the corresponding 
properties of the Lie bracket. The Leibnitz rule follows from the Leibnitz rule for differentials. 

The D.b. (|^T|) preserves the class of gauge-invariant functions. Suppose f,g are gauge- 
invariant. Then df,dg G SVect(M, /i). Recall (Sec. 12. 3j) that if / and g are gauge-invariant, 
then so are their differentials d/* and dg' . It follows that [df, dgf is also gauge-invariant. Now, 
under a gauge transformation A' = A + dA , 



{f,g}{A') = {A',[df,dg]{A')) = {A',[df,dg]{A))= J A',{[df,dg]{A)yp 

{f,g}{A') - {f,g}{A) =-J Adim,dg]{A)yp)d^x = (32) 

Thus H31() is a gauge invariant Poisson bracket. 

A gauge-invariant observable f{A) defines canonical transformations on the Poisson manifold 
Q* via the p.b. Suppose g{A) is any observable, then its Lie derivative under the flow generated 
by / is Cvf9{A) = {f,g}iA) = {A, [df,dg]) (see Sec-B- 

Example: The p.b. of two moment maps PuiA) and PviA) is 

{pu,Pv}iA) = P[u,v]iA) = J^Ai{u^djv' - v^dju')p (33) 
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If u and V are volume preserving, then so is [u, v\ ; therefore, if and are gauge-invariant 
functions of A, then so is ^J^■[u^v]■ Moreover, the canonical transformation generated by the 
moment map ^i^ is just the Lie algebra coadjoint action (see ^ and Sec. |3J) 

C-v,J[A) = {fiuJ}{A); Cv,^A = a<A. (34) 



{A{x),A,{y)} = 6\x-y) 



(35) 



3.2 Poisson Brackets of Gauge Fields 

Equivalence classes of gauge fields are coordinates on our Poisson manifold G* = SVect(M, /i)* . 
Thus, an explicit formula for the 'fundamental' p.b. between components of the gauge field 
Ai{x) is useful. This will also facilitate a comparison with electrodynamics. We will show that 
the p.b. between gauge field components is 

My)^p'Hy) - A^{x)^p~\x) 

where the derivatives act on everything to their right. Though this formula for {Ai{x), Aj(y)} 
looks a bit complicated, the rhs is linear in gauge fields. Thus, our Poisson algebra is actually 
a Lie algebra like the Lie algebra of angular momenta {Li, Lj} = CijkLk . 

Recall that electrodynamics is based on the Heisenberg algebra between the spatial compo- 
nents of the gauge field and the spatial components of the electric field. In two spatial dimensions 
we have (before any gauge fixing) 

{M^,t),E^{y,t)} = 6^{^-y)d{ 
{Aiix,t),Ajiy,t)} = 

{E'i^,t),EHy,t)} = (36) 

While in electrodynamics the electric field is canonically conjugate to the gauge field, this is 
not the case in our theory. The components of the gauge field in our theory obey p.b. rela- 
tions with each other, without being canonically conjugate. This is not unusual. For instance, 
the components of angular momentum form a closed Poisson algebra though none of them is 
canonically conjugate to another. This is a generic feature of degenerate Poisson manifolds 
where canonically conjugate variables can only be chosen on individual symplectic leaves (the 
concentric spheres in the case of angular momenta). 

Gauge fields obeying p.b. not involving the electric field are not alien to conventional Yang- 
Mills theory. For example, in a coordinate system where initial values of fields are specified on a 
null cone at past time-like infinity, the transverse components of gauge fields satisfy p.b. among 
themselves as shown by Rajeev and Turgut[2linj 

1 

where are transverse angular coordinates, qij is the round metric on 5^, ii is a radial 
coordinate and a, b are color indices. In fact, our original motivation for studying the dynamical 
system in this paper was to find a toy-model that shared this feature with Yang-Mills theory. 

Now we will establish (|35j) using (|31j) and the relation between the p.b. of functions and 
those between the 'coordinates' Ai{x) 

{MiA) = /^^-^^.{^.(-),^.(.)}^^ (38) 



{Al{z,R),A'j,{z',R')} = -5-,5lq,,{z)5{z - z') sgn {R- R') (37) 
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We rewrite {f,g}{A) from ()31() to make it look like this 





{- 


^9 ^ 


(- )] 


Ml'' 


\p5A,) 


-SA^ 


\pSAj\ 



d^x(fy5^{x-y)Aj{y) 
d^x(fy6^{x-y)Aj{y) 



6f d / I 6g 



dP'xcPy5'^{x — y) 
Finally we read off the desired expression 1)35^ . 



5Ai{x) dy^\p{y)6Aj{y) 
9 1 Sf 6g 
dtipiy) SA{x) 6Aj{y) " 



f ^ 9 



My)Q-iP''{y)-Mx)-Q-P 



6Aiix) 6Aj{y)' 



(39) 



4 Structure of Poisson Algebra of Observables 



In this section we give the canonical action of SDiff (M,;u) (generated via p.b.) on the Poisson 
manifold G* = SVect(M, /x)* (see also^ Ref. ^S]). Symplectic leaves are coadjoint orbits of 
SDiff(M, ^). The Poisson algebra of gauge invariant observables is degenerate. The enstrophy 
invariants In of hydrodynamics are an infinite number of Casimirs. They are constant on the 
coadjoint orbits. Then we try to identify some of the simpler coadjoint orbits. We show that 
there is always at least one single-point orbit, that of the pure gauge configuration. If M has 
non- vanishing first cohomology, then we show that there are finite dimensional symplectic leaves 
lying inside H^{M) \ {0}. This also shows that /„ could not be a complete set of coadjoint 
orbit invariants. If M is simply connected, we show that the only single-point orbits consist of 
the configurations for which dA/fi is constant. For simply connected M , we also argue that 
the orbit of [A] € G* for which dA/p, is not constant, is infinite dimensional. We identify the 
isotropy sub- algebra and tangent space of such an orbit and give an example where dA/fi is 
circularly symmetric. However, this analysis is far from complete. It would be useful to find a 
nice coordinate system on these orbits and get the symplectic structure in explicit form so as to 
study hamiltonian reduction. 



4.1 Coadjoint Action of SDiff (M,/i) on Poisson Algebra 

The Poisson manifold G* = n^{M)/d^f{M) = {A e n^{M) \A ~ A + dA} carries the coadjoint 
action of SDiff(M, /i). For example, on a simply connected region M, G* may be identified 
with the space of scalar functions / = (dA/fi) = (B/p) on M. The coadjoint action is the pull- 
back action of volume preserving diffeomorphisms 4> G SDiff (M,//) on functions (p*f = / o 0. 
The action on equivalence classes of gauge fields is also the pull back = For 

infinitesimal (f){t) = 1 + ut we get the Lie algebra coadjoint action of n G SVect(M, ^) on G* 

adlA = -CuA (40) 

This action of SDiff (M, p) on G* is canonical i.e. there is a function on G* (a gauge-invariant 
observable) that generates the coadjoint action via the p.b. The generating function is the 
moment map pu{A). To see this, suppose u is a volume preserving vector field, then 

{pu{A),Aj} = -{CuA)j. (41) 
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To show this we begin with the rhs of (|41() 




(42) 




y) 



= -p ^di{Akpu') - Aidku' = -u'diAk - Aidku' = -{CuA)k 



(43) 



We integrated by parts (yl = on dM) and used di{pu^) = 0. Thus, we obtain (|1T|) . 
4.2 Center of Poisson Algebra 

From Sec. 14.11 the symplectic leaves of SVect(M, /i)* are homogeneous symplectic manifolds, 
identified with coadjoint orbits of SDiff (M, p) . Is the Poisson algebra of functions on SVect(M, p)* 
degenerate? What are the Casimirs that constitute its center? Casimirs are unchanged under 
canonical transformations. So they are constant on symplectic leaves. A complete set of such 
Casimirs would allow us to distinguish between distinct leaves. Since the symplectic leaves are 
coadjoint orbits of ST)iS{M, p) , it suffices to find observables that commute with the moment 
maps pu{A) , which generate the coadjoint action. 

We first observe that the center of the Poisson subalgebra of linear observables PuiA) is 
trivial, {pu, Pv}iA) vanishes for all A iff [u,v] = 0. But the center of SVect(M, /x) is trivial. 
This is seen by writing u and v in terms of their stream functions ipuji^v Taking some simple 
choices for in the condition [u,v] = will imply that ip^ = 0. Thus none of the Pu{A) lie in 
the center of the Poisson algebra of gauge- invariant functions, since they do not even commute 
with each other. Thus we need to look elsewhere to find the Casimirs of our Poisson algebra. It 
turns out that the charges /„ = J{dA/p)^p are central observables. The monomials [dA/p)^ 
can be replaced by any scalar function of (dA/p). We found that In are Casimirs by showing 
that they Poisson commute with the moment maps. However, this involves lengthy calculations. 
For eg., in El we show that commute with each other and in Owe show that I2 commutes 
with Pu{A) for uniform p. 

However, we found that In are closely related to the enstrophy invariants of hydrodynamics, 
which are known to be conserved quantities in Eulerian hydrodynamics. There is a simple 
argument in Ref. jin] based on earlier work|2Tl 1201 122] that leads to the conclusion that In 
are constant on coadjoint orbits. The argument is that the action of SDiff (M,/x) is merely to 
change coordinates in the integral defining /„ . Since this integral is independent of the choice of 
coordinates. In must be invariant under the coadjoint action. Though I„ for n = 1,2,3,... are 
all Casimirs, they are not a complete set of orbit invariants. Their level sets do not necessarily 
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distinguish the coadjoint orbits (see Sec. 14.3(1 . For some remarks on additional invariants, see 
Sec. 9 of Ref. JJ.^. Our explicit calculations of {In, fJ'u} given in El and [0 suggested to us how 
/„ could be modified in order to get an independent infinite sequence of conserved quantities 
(not Casimir invariants) for our choice of hamiltonian; see Sec. 16.31 

4.3 Finite Dimensional Symplectic Leaves in H^{M) 

The simplest symplectic leaf in the dual of the Lie algebra of the rotation group is the origin 
of angular momentum space {Li = 0}. It consists of a single point. Of course, one can also 
characterize the orbit {Li = 0} as the zero set of the Casimir = 0. Can we get a similar 
explicit characterization of the simplest symplectic leaves of SVect(M, /x)* ? 

Orbits of Closed One Forms and Zero Set of Casimirs In 

The simplest symplectic leaf in SVect(-/Vf, /i)* should be the orbit of exact 1-forms. Suppose 
A = dA is an exact 1-form. Then under the coadjoint action of n G SVect(M, /i) it goes to 
A' = A + adlA 

A'i = Ai- iCuA)i = diA - {djA){diU^) + u^djdiA = di{A - u^d.A) (44) 

We see that a pure gauge is mapped to a pure gauge under the Lie algebra coadjoint action. Thus 
the tangent space to the orbit of exact 1-forms is trivial. So pure gauges form a single-point 
orbit. 

After the pure gauges, the next simplest configurations are closed but inexact 1-forms dA = 
A ^ dA. Suppose M is a manifold with non- vanishing first cohomology. What is the orbit 
of an element of H^{M)1 In particular, is the orbit a finite dimensional manifold? Does the 
orbit lie within H^{M)1 The answers to both these questions is affirmative. Suppose A is an 
exact 1-form, F = dA = 0. Under the Lie algebra coadjoint action. A' = A + ad^A = A — CuA. 

dA' = d{A - CuA) = -d{diu + iud)A = -diudA = (45) 

The tangent space to the orbit of a closed 1-form contains only closed 1-forms. Since the space 
of closed 1-forms on a two-dimensional manifold is finite dimensional, the tangent space to the 
orbit must be finite dimensional. Moreover, we have shown that the exact 1-forms form a single 
point orbit by themselves. Thus, at the infinitesimal level, the orbit of a non-trivial element of 
H^{M) must lie within H^{M) \ {0}. This view will be reinforced by considering the zero set 
of Casimirs 

In = I {dA/fir^i = (46) 

/„ are constant on coadjoint orbits. Thus, the orbits must be contained in their level sets. 
I2n is the integral of a positive quantity and therefore vanishes iff dA = 0. If dA = then 
l2n+i is also zero. Thus the level set /„ = is the space of closed 1-forms on M . Thus, the 
coadjoint orbit of a closed one form must lie in H^{M). Since we already established that the 
pure gauges form a single-point orbit, this means In cannot be a complete set of Casimirs if 
H^(M) is non-trivial. 
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For example, if M is the plane, then iif^(R^) consists only of the equivalence class of pure 
gauge configurations and the zero set of contains only one single-point orbit [^4] =0. The 
same is true of the 2-sphere which has trivial first cohomology. For the 2-torus i7^(T^) ~ . 
For , there must be symplectic leaves which are submanifolds of \ {0} . It is interesting to 
find the orbits of cohomologically non-trivial gauge fields more explicitly as well as the induced 
symplectic structure. However, we do not investigate this further since the Hamiltonian we pick 
H66|) vanishes on closed 1-forms. There is no interesting dynamics on the finite dimensional 
symplectic leaves we have described above. Therefore, we turn to the case where M is simply 
connected and try to characterize the orbits of gauge fields that are not closed. 

4.4 Symplectic Leaves when M is Simply Connected 

In the case of the rigid body, symplectic leaves other than Li = are concentric spheres of non- 
zero radius, all two dimensional symplectic manifolds. These leaves may also be characterized as 
non-zero level sets of . By analogy, we would like to find the orbits in SVect(M, ^)* of 1-forms 
that are not closed. They must lie within non-zero level sets of . Can we say something more 
about them, such as whether they are finite dimensional? We address these questions below, 
assuming M is simply connected. 

By Poincare's lemma, all the information in a gauge field on a simply connected M can be 
stored in the 2-form field strength, F = dA. Then, the dual of the Lie algebra SVect(M, /x)* 
may be identified with the space of scalar functions on M 

dA 

SVect(M,/x)* = {/ = —} (47) 

The coadjoint action of (/> G SDiff (M, /i) is just the pull back 

Ady = <p*f = fo(t> (48) 

and the Lie algebra coadjoint action of u G SVect(M, ^u) is 

/ ^ / + adlf = f-Cuf = f- u'dif (49) 

Thus the coadjoint orbit of / and the tangent space to the orbit at / are 

Of = {/o</. j E SDiff (M,/i)} 
TfO = {/ - u'dif I u G SVect(M,/i)} (50) 

The picture that emerges from our analysis below, is that there are two types of orbits when 
M is simply connected. Orbits of the first type contain only a single point, namely a constant 
function / = dA/fi = c. In the case of the plane with the uniform measure, the only one-point 
orbit is {/ = (dA/fi) = 0} , consisting of the pure gauge configuration. The orbits of the second 
type are all infinite dimensional. They are the orbits of non-constant functions /. In the case 
of the plane with uniform measure, we expect an infinite number of such orbits of the second 
type, each contained within a level set of {/n,"- = 1,2,3, .. .}, though we do not rule out the 
existence of more than one such orbit in any one level set of the /„ . 
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Single Point Orbits 



We show that constant functions / = (dA/fi) = c are the only single point coadjoint orbits 
in the dual of the Lie algebra SVcct(M, /v,)* when M is simply connected. They lie within the 
level sets /„ = Vol(M,/x) c'' where Vol(M,^) = J^l^- 

Suppose / is a constant. Then Cuf = 0. So the constant functions form single point orbits. 
The case / = corresponds to the closed and exact 1-forms A which we already identified as 
a single point leaf if M is simply connected. Non-zero constant functions are not admissible 
elements of SVect(M, jj,)* if M = and fi is the uniform measure. But they are allowed if M 
is compact or if ^ at infinity on non-compact M . Note that In = J f^n = Vol(M, /i) c" if 
f = c. So the constant functions lie in the level sets with /„ = Vol(M, /i) c". 

We can go one step further and show that constant functions are the only single point orbits 
if M is simply connected. Suppose / is a single point orbit. Then must vanish for all 
volume preserving vector fields u. Since M is simply connected, any such vector field can be 
written in terms of a stream function = p^^e^^djip. Then denoting derivatives by subscripts, 

Cuf = ^ fxipy = ^xfy (51) 

This must be true for all stream functions ip . Taking ijj = x and = y successively tells us 
that / must be independent of both x and y and hence a constant. We conclude that the only 
single point orbits in Q* = SVect(M, /i)* are the constant functions f = (dA/p) = c when M 
is simply connected. 

Orbit and Stabilizer of Non-constant Element of SVect(M, /li)* . 

The stabilizer algebra Stab(/) or isotropy subalgebra of a function / is the set of all volume 
preserving vector fields which leave it fixed under the coadjoint action. The cosct space 
SVect(M, /n)/Stab(/) then has the same dimension as the tangent space to the orbit containing 
/. For u to be in Stab(/) we need £„/ = 0. Since M is simply connected we can express u in 
terms of its stream function = p~^e^^djX- The condition = becomes 

e'-'{d,f){djx) = ^ Xxfy-Xyfx = (52) 

which in vector notation^ says that V/ x Vx = 0. Colloquially, the gradient of x must be 
everywhere parallel to the gradient of /. x = c/ is clearly a solution for any real number c, so 
the isotropy subalgebra is at least one dimensional, provided / is not a constant. The product 
of two solutions as well as real linear combinations of solutions are again solutions to this linear 
PDE. To better understand the general solution of this PDE, let us first consider the specific 
example of a function f(r,9) that is circularly symmetric. 

Example: Stabilizer and Orbit of Circularly Symmetric Function 

For example, let / = (dA/p) be a non-constant function depending only on the radial coordinate 
r = yx^ + y^ • Let us call its orbit by the name O / . The gradient V/ points radially, so to 

We emphasize that this equation is independent of any metric on M . 
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speak. The isotropy subalgebra of any sueh funetion / is the space of stream functions % with 
V/ X Vx = 0. The solutions are stream functions xi'>') tti^t are independent of 9 

Stah{f) = {x{r,e)\d0X = O} (53) 
In this case, the isotropy subalgebra is infinite dimensional. The coset space 

SVect(M,^)/Stab(/) = {^(r, ^)|V' ~ V + x(^)} (54) 

has the same dimension as the tangent space to the orbit TfOf at /. Though the stabilizer is 
infinite dimensional, the orbit of /(r) is infinite dimensional as well. For eg. the coset space 
can be parameterized using an infinite number Fourier coefficients 

oo oo 

^(r, e) = J2 i^n^ (r) COS nd+J2 '^n^ (r) sin ne (55) 

n=l n=l 

Wc omitted the 6 independent additive term in order to quotient out by Stab(/) . 

To summarize, the isotropy subalgebra of a non-constant circularly symmetric function con- 
sists of all stream functions that are circularly symmetric. Moreover, the tangent space to the 
orbit may be identified with the coset space in which two stream functions are identified if they 
differ by one depending on r alone. 



Non-constant Functions have Infinite Dimensional Orbits 

Using the circularly symmetric example as a guide, we can characterize the isotropy subalgebra 
of a general non-constant function. The condition for x to be in the isotropy subalgebra is 
fxXy ~ fyXx = . We will describe the set of all solutions x to this linear PDE in a region of M 
where df is never the zero 2-form. Such a region is guaranteed to exist by an analogue of the 
inverse function theorem, since we assumed that / was not a constant. Our answer is that the 
isotropy subalgebra of / consists of all stream functions x which are constant along level sets 
of /. The level sets of / are necessarily one dimensional curves in such a region. 

To see this, we argue as follows. In a region where f{x,y) has no local extrema, the level 
sets of / are one dimensional. These level curves foliate the region. Pick a coordinate 0{x,y) 
along the level curves. Then def = /e = 0. Also pick a coordinate R 'transversal' to the level 
curves of /. What this means is that Or and Oq should be linearly independent at each point 
{R, Q) . There is a lot of arbitrariness in the choice of these coordinates, and one certainly does 
not require any metric to define them. The volume form /x is non-degenerate, so the volume 
enclosed by the mini-parallelogram determined by the tangent vectors Or and Oq is non-zero. 
Denoting partial derivatives by subscripts, we get 

Hidu, de) = pdx A dy{dR, do) 

= -p{dx <Sidy-dy<Si dx){—dx + -^dy, —d^ + -^dy) 

1 1_\ ^ / RyOj-RxSy 

2 \RxQy Ry&xJ 2 \ RxQyRyQx 

Since p is non-vanishing, this means 

RyQx - Rx&y + (57) 




16 



The condition for x to lie in the isotropy subalgebra becomes (using /e = 0) 

fxXy-fyXx = UrRx + fe&x){XRRy + XeQy) - X <^ y 

= fRRxiXRRy + Xe0j/) - fRRyiXRRx + XS^x) 

= fRXeiRxOy - RyOx) (58) 

We have already shown that RxQy — Ry&x / 0. Moreover, / is non-constant and /e = so 
we must have /r 7^ 0. We conclude that x& = 0. In other words, the stream functions x 
the isotropy subalgebra are constant on level sets of /. 

Stab(/) = { Stream functions x constant on level sets of /} (59) 

It is clear that the space of such stream functions is closed under linear combinations and 
also under multiplication, as we observed earlier. The tangent space to the orbit at / is then 
identified with the coset space 

TfOf = SVect(M, /x)/Stab(/) = {7p\ip ~ V' + X, X constant on level sets of /} (60) 

Both the stabilizer and the orbit are infinite dimensional, provided / is not constant, for the 
same reason as given in the example where / was circularly symmetric. Thus, non- constant 
f = dA/fj, have infinite dimensional orbits. 



5 Hamiltonian 



Having specified observables (gauge- invariant functions), their p.b. and phase space (coadjoint 
orbit of SDiff (M, /i) ) , we need to pick a gauge invariant hamiltonian. 



5.1 Hamiltonian Leading to Ideal Hydrodynamics 

There is more than one interesting way to pick a hamiltonian. The classic choice leading to ideal 
hydrodynamics, requires a positive metric gij on M which defines a positive, symmetric, inner 
product on ^ = SVect(M, /i) 

{u,v)g= / g{u,v) fi= / giju'v^ n. (61) 
Jm Jm 

The hamiltonian of Eulerian hydrodynamics is then 

1 1 f ■ ■ 

H{u) = -{u,u)g = - J gij p d^x (62) 

There is no a priori relation between the metric g and the volume form fi. Qg = ^/det gdx^ A dx'^ 
need not equal fi. If /x = $7^, the theory is particularly natural as well as non-trivial. The 
hamiltonian defines an inertia operator / (generalization of the inertia tensor of a rigid body, 
seeEl Ref. [13) from the Lie algebra G = SVect(M,/x) to its dual G* = (M) / dft^ (M) . 
Suppose u, V are volume preserving vector fields. Then the inertia operator is defined by the 
equation {Iu,v) = {u,v)g. In other words. 



{Iu)jV^fi = / gtjU^v^ IX =^ I {{Iu)j - giju') p = (63) 
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Since this is true for an arbitrary volume preserving v, {{Iu)j —giju^) must be a total derivative: 
{Iu)j — QijU^ = djA for some scalar function A which vanishes on the boundary of M or at 
infinity. We see that the metric does not determine the 1-form lu uniquely, but rather up to an 
exact 1-form. Thus, the image lu of a volume preserving vector field u £ G is an equivalence 
class in G* = gauge fields modulo gauge transformations. A coset representative is given by the 
1-form Ai = QijU^ . The equation of motion is the well-known Euler equation 

= -Vuu -Vp; CufJ. = (64) 

Here V^w is the covariant derivative (with respect to gij ) of u along itself. The two equations 
fix the pressure p{x, t) up to an additive constant. The Euler equations of hydrodynamics have 
a geometric interpretation^]. The configuration space of a volume preserving fluid flowing on 
the manifold M is the group of volume preserving diffeomorphisms SDiff(M, ^u). The inner 
product (|61|) on the Lie algebra of this group can be extended to a right-invariant metric on the 
whole group by right translations by group elements. Then, by the least action principle, the 
time evolution of the fluid is given by geodesies on SDiff(M, //) with respect to this metric. 



5.2 Hamiltonian as Magnetic Energy 

Now we propose a hamiltonian different from that of Eulerian hydrodynamics. It is a magnetic 
energy, natural from the point of view of Yang- Mills theory. The result will still be a theory of 
geodesies on SDiff (M, /i) , but with respect to a different right-invariant metric on this group 
than that implied by ideal hydrodynamics. Suppose M is a 2-dimensional surface with volume 
form fi and metric g^j . The metric did not play any role so far since the phase space and 
poisson structure are independent of it. But to specify the dynamics, we need the metric. Any 
two dimensional manifold is conformally flat, so the information in g^j is encoded in its volume 
form = ^/gdx^ A dx"^ , where g = det gij . We do not assume that /j, is equal to i7g . Indeed, 
if that is the case, the dynamics is trivial. By analogy with Yang-Mills theory, we pick the 
manifestly gauge-invariant magnetic energy 

as our hamiltonian. It can be written in a variety of equivalent ways 

We find the last formula most useful in calculations. Here a is the positive scalar function on 
M given by the square of the quotient of the two volume forms. 

a = ifi/ngf = p^g (67) 

We will see that for the dynamics to be non-trivial, a must not be constant. The hamiltonian 
in (|66|) is the simplest choice that is gauge-invariant, quadratic in gauge fields and involves two 
derivatives. It is easy to see the equivalence of the various formulae for the magnetic energy in 
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Let e'^ = g'^g^^eu, e'^ = = 1/5- Moreover, ge'^ = e'K The Hodge dual of F is 

the scalar 



* F = ^{dx^ A dx^)-Fii = -Jge^^Fij = V Fj.ej, - — 



F A *F is the 2-form [B"^ / ^)dx^ A dx^. ((F A *F)/r2g) = B"^ / g is a scalar. Hence 

(F A *F)/J^, = (*F)2 = (F/f7,)2 = = = (F/;u)V (69) 

Moreover, Fi^F'^ = 2{Fi2F^^). But F^^ = ^i^^^^F^, = {g^^ g^^ - g^^ g^^)Fi2 = g-^Fi2 = B/g. 
Thus FijF"^^ = 2B'^/g. This shows the equivalence of all the forms of the hamiltonian given in 

The hamiltonian ()66() involves only the spatial components of the field strength since there is 
no time component for the gauge field in our theory. So there is no analogue of electric energy. 
Moreover, our theory is not relativistically covariant unlike electrodynamics. The dynamics 
determined by ()66() is not equivalent to Eulerian hydrodynamics, though both theories share the 
same phase space. To see this, note that ii Qg = fj,, then ()66() reduces to the Casimir I2 and 
has trivial dynamics while the hamiltonian H62|) continues to have non-trivial dynamics. 

The differential of the hamiltonian is ( qa^x) ~ ~£^''dj5'^{x — y)) 

{dH{A)y = = -e^^d,{Ba/p), (70) 

For each A, dH{A) is volume preserving di{pdH^) = e^^ didj{Ba/ p) = 0. H generates a 1- 
parameter (time) family of diffeomorphisms of a coadjoint orbit of SDiff(M, ^u), which serves 
as the phase space. The hamiltonian vector field Vh is a vector field on a coadjoint orbit O. 
At each tangent space Tj^^^,)]© to an orbit, Vh([^]) is given by the coadjoint action of the 
Lie algebra element dH{A) , Vh = o.d'^H(A) ■ "^^^ change in an observable / under such an 
infinitesimal canonical transformation is the Lie derivative with respect to Vh 



df{A) 



-Cv„f = ad*^Hf = {H, /} = (A, [dH, df])= j Ai[dH, dff p (71) 



dt 

Explicit equations of motion in local coordinates are given in Sec. 



5.3 Inertia Operator and Inner Product on Q* from Hamiltonian 

Recall (Sec. I5.1|l that the hamiltonian of ideal hydrodynamics (|62|) defines an inner prod- 
uct on the Lie algebra of volume preserving vector fields via the inertia operator. Here, 
the magnetic energy defines a positive inner product {■,-)g* on the dual of the Lie algebra, 
Q* = (M) / dOP {M) , via an inverse inertia operator h, obtained below. This inner prod- 
uct is degenerate if M has non- vanishing first cohomology. Away from degeneracies, one can 
in principle invert h to obtain an inertia operator and an inner product {■, ■)g on the Lie al- 
gebra G = SVect(M, /i). It should be possible to extended to a metric on the group 
G = SDiff (M, p) by means of right translations. The equations of motion in Q* = SVect(M, p)* 
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determined by H should be the projection of the geodesic flow on G = SDiff(M, //) with re- 
spect to this right-invariant metric. In this sense, the dynamics determined by the magnetic 
energy describes geodesies on SDiff(M, /x). Eulerian hydrodynamics also describes geodesies 
on SDifl'(M, yu), but with respect to a difl'erent right-invariant metric lS;- Here, we obtain the 
inverse of the inertia operator and the inner product on SVect(M, ^)* explicitly. Let us begin 
by writing the hamiltonian as a quadratic form on gauge fields. The result is 

H = - Aih'^ Ai p (fx with : 



P 



{dj-)dk + -djdk 
P P . 



(72) 



h^^ plays the role of inverse inertia operator, mapping equivalence classes of gauge fields to 
volume preserving vector fields by raising the index Aj i-^ h^'' = idH{A)y . To get (|72)) . note 
that H may be written in terms of its differential, 

H = ^iA, dH{A)) = \pdH[A) = \j A, {dH{A)y p d^x. (73) 

This is easily checked by integrating by parts assuming S ^ on dM: 



_ 1 
2 

As a consequence. 



H = \ j {B\/p)d^x = ^ 1 {Ba/p)d,A,d^x = '-^ J d,{Ba/p)Ad^ 

AidWpd^x (74) 



H = \ I AiihAYp. where {hA)' = h'^ Aj = dW = —dj{Ba/p) (75) 



Writing this out we get. 



(^hAy = L^d.(^o^Ar 



^ ^idjia/p))dk + ia/p)djdk^Ai (76) 

From this we read off the inverse of the inertia operator (|72j) . Through h, the magnetic energy 
naturally defines a symmetric positive inner product on the dual of the Lie algebra, G* = 
(M) / dQ^ (M) . This inner product may be written in the following equivalent ways 

{A,A)g* = (A, hA) = Pj^A^A) = j Ai^Aj pd^x = j B B {^d^x (77) 

Here B,B are the magnetic fields corresponding to A, A. Positivity of the inner product is 
ensured since (j{x) = (p/Qg)'^ = {p^ / g) > and the integrand is positive 

{A,A)g, = 2H{A) = J {B'^a/p)d^x > (78) 

Symmetry of the inner product is shown by establishing the last equality in (|77)) : 

j Ai{hAypd^x = e'^ j Aidj{Ba/p)d^x = -e'^ j djAi{Ba/p)d^x 

\BBa/p)d^x. (79) 



Inner product (|77|) is degenerate precisely on gauge fields in the first cohomology of M, H^{M) = 
{B = 0,^7^ dK} . If M is simply connected, H77|) is a non-degenerate inner product on 
Q* = SVect(M, //)* . Inverting it gives an inner product on Q which may be extended by right 
translation to a right invariant metric on SDiff (M, p) . 
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6 Equations of Motion 



Before working out the equations of motion in detail, it would be prudent to convince ourselves 
that unlike the Casimir I2 , the magnetic energy H does not lie in the center of the Poisson 
algebra for non-constant a . We show this in ^ by exhibiting a moment map that has non- 
vanishing p.b. with H . 



6.1 Time Evolution of Magnetic Field 

The time evolution of the magnetic field is given by the equation 

B = {H,B} = e'^di{B/p)dj{Ba/p) = V{B/p) x V{Ba/p) 



(80) 



where a = p^ / g ■ It is manifestly gauge-invariant. Though the hamiltonian is quadratic, (|8()|) 
is non-linear since the p.b. of gauge fields ()35() is linear in gauge fields. The 'interactions' are 
partly encoded in the Poisson structure and partly in the hamiltonian, so to speak. (jHOJ is an 
analogue of the Euler equation L = L x $7, L = IQ for the rigid body, (see^. Indeed, the 
equation of motion for the magnetic field can be regarded as the Euler equation for the group 
SDiff (M, p) with respect to the metric ()77() defined by the hamiltonian ()66|). The time evolution 
of the magnetic field can also be written as an equation of motion for the field strength F = dA 
defined in (EH) 



F 



1 



e'^di{F/p)dj{Fa/p) eudx^ A dx^ 



(81) 



The formula 1)80(1 for B is obtained using the p.b. formula (|31|) using the relation sAj{y) 



^jk dS [z^ y) ^ g^^^yumiig yl = on dM . After some integration by parts. 



dB 
'dt 



{H,B} 



p-\diAj) + Ai{djp-^) + -Aidj 



diiBa/p) 



(82) 



Expanding out derivatives of products, eliminating A in favor B = e^WiAj and after a lot of 
cancelations, one arrives at the advertised evolution equation. 



6.2 Time Evolution of Moment Maps 

We have already argued that the Casimirs I„ = j{dA/ p)"^ p are constant on symplectic leaves. 
Thus, they are conserved under time evolution = {H, /„} = 0, independent of the choice of 
hamiltonian. As for the moment maps pu , we can show using ()31() that 



{H,pu} = / d'^xA 



dpujA) 
dt 

For p = 1 this reduces to 

dpu{A) 
dt 



dk{Ba/p) 



dk{Ba) 



(83) 



(84) 



The rhs is gauge-invariant even though the gauge field appears explicitly. 
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6.3 Infinitely Many Conserved Charges in Involution 



We find an infinite set of conserved charges 

dHr 



dt 



{H,Hn} = 0; n = 1,2,3,... 



(85) 



for a uniform fi {p = 1) and gij an arbitrary metric^. We suspect that a similar result holds 
for non-uniform fi. The conserved quantities are 



(86) 



where a = (fi/flg)'^ = (l/g) and g = det gij. Furthermore, we find that Hn are in involution 
{Hm,Hn} = 0. The hamiltonian is H = \H2. The presence of infinitely many conserved 
charges in involution suggests an integrable structure underlying the dynamics determined by 
the magnetic energy. The hamiltonian of ideal hydrodynamics (|f)2|) is not known nor expected 
to have any conserved quantities besides the /„ and their close relatives, which are constant 
on coadjoint orbits. What is remarkable about Hn is that they are not constant on coadjoint 
orbits but still conserved quantities for the magnetic energy hamiltonian ()66() . To establish this 
we investigate the time evolution of Hn ■ 



dHn 
dt 



{H,Hn} = J Ai 



6H 



dj p 



-1 



5H„ 
5A 



SH 



SA, 



-dj p 



^_,5HX 
6aJ 



d^x 



(87) 



Using 



We get 



5Hn 



ne'^dj[(j{B/p) 



n-l 



dHn 
"dt 



For p = 1 this becomes 



ne- 



^^e'^ / Ai 



dk{aB/p)d,{p-^di{a{B/pT-^)) 



-dk{a{B/pT~^)d,{p-^di{aB/p)) 



d^x 



dHn 



A, 



dk{aB)djdi{aB 



dkiaB^-')djdiiaB) 



d^x 



(90) 



In the gauge = , A2 = A, B = dxA this may be written as 



dHn 
dt 



n J dxdyA 



dx{Ba)dxy{B''''a) - d,{B^~'a)d,y{Ba 



dy{Ba)di{B''''a) + dy{B'^-'(j)di{Ba) 



(91) 



We prove in^that Hn, n = 1, 2, 3, . . . are conserved by showing that the rhs of H91|) vanishes. 
Having found an infinite number of conserved quantities we wanted to know whether their p.b. 



*o" = (? 1 9 must not grow too fast at 00 and the magnetic field must vanish at 00. 
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generates new conserved quantities. In our case, we discovered (again for p = 1), that the 
conserved quantities Hn are in involution, i.e. they mutually Poisson commute {Hm,Hn} = 0. 
The proof of this is somewhat lengthy and is relegated to El 

Intuitively, Hn are independent of each other since they are like average values of different 
powers of B . Furthermore, Hn are independent of the Casimirs Im ■ Im contain no information 
about the metric gij while H^ depend on the metric via a . For example, we show in IdI that 
H2 is not a Casimir. More generally, it would be nice to prove that I„ and H^a are functionally 
independent by showing that on every tangent space to an orbit, the cross product of their 
gradients is non-vanishing. 

7 Static Solutions 

7.1 Zero Energy Configurations 

The magnetic energy is H = Jj^j{dA/ n)'^afi = J(B'^ / p)ad?x where a = (/x/il^)^ = {p^ / g) > 0. 
Thus, the energy is non-negative, H >0. Moreover, if A is closed, the energy automatically is 
a global minimum dA = ^ H = 0. Moreover, since H is the integral of the square of dA, 
weighted by a positive function, closed gauge fields are the only configurations with zero energy. 
Any such closed gauge field is a static solution to the equations of motion (irrespective of p and 

9ij) 

%\dA=o = {HJ}\dA=o = (92) 
at 

As discussed in Sec. 14.31 the closed gauge fields constitute one or more (according as M is simply 
connected or not) symplectic leaves of the Poisson manifold. They correspond to the zero set 
of Casimirs In = JidA/py^p = 0. Thus, the hamiltonian vanishes on all the symplectic leaves 
with In = and there is no interesting dynamics to speak of. If M is the plane, then the pure 
gauges are the only ones with zero energy. The leaf /leaves with /„ = are the analogue of the 
= point at the origin Lj = of the angular momentum Poisson manifold. At that point, 
the energy E = J^i^if^^i of the rigid body vanishes as well. 

It is interesting to find minima of energy on more interesting symplectic leaves. Initial 
conditions determine which symplectic leaf is the phase space of the theory. The general problem 
of finding the minimum of energy on a given symplectic leaf (perhaps specified through values 
of invariants such as /„) is potentially quite interesting and difficult. One would first have to 
find which gauge fields or magnetic field configurations satisfy the constraints and lie on the 
specified leaf. This is similar to the problem we solved in the large N limit of 2d QCD where we 
found the minimum of energy on the symplectic leaf with baryon number equal to one jlll ITU] . 
Here we do the opposite, find a few static solutions and then determine which orbit they lie on. 

7.2 Circularly Symmetric Static Solutions 

Suppose M = . If both the volume forms p and i7g are circularly symmetric, then any 
circularly symmetric initial magnetic field will remain unchanged with time. To see this, note 
that a = p^ / g depends only on the radial coordinate r. Suppose that at t = 0, B{r,t = 0) 
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depends only on r. The initial value problem for B given in H8(J|) is 

B = V{B/p) xV{Ba/p) (93) 

Due to circular symmetry, both the gradients point radially and their cross product vanishes. 
Thus B = and we have a static solution B{r) ! 

We already met the B{r) = static solution before, it lies on a one point symplectic leaf 
of pure gauge configurations. However, the static solutions corresponding to non-constant B{r) 
lie on infinite dimensional symplectic leaves which we found previously (Sec. 14.41) . The values of 
Casimirs on the orbit containing a circularly symmetric static solution B{r) are 

In = 27r {B{r)/p{r)Y p{r) r dr (94) 
Jo 

By a judicious choice^ of B{r) one should be able to find a static solution B{r) for given p{r) 
that lies on an orbit with practically any desired value for the invariants I„ . More generally, by 
an argument similar to the one given in Sec. 14. 4( we see that magnetic fields for which (B/p) 
and (Ba/p) have common 1-dimensional level sets, are static solutions of 



7.3 Some Other Local Extrema of Energy 

For a gauge field to be an extremum of energy on a given symplectic leaf, the variation of energy 
in directions tangential to the leaf must vanish. There is no need for variations in other, let alone 
all, directions to vanish. However, though it is not necessary, if [A] is such that all variations of 
the hamiltonian H{A) vanish, then, A must be a local extremum of energy. Such an extremum 
= is automatically a static solution to the equations of motion for any gauge-invariant 
observable / 

I = {HJ} = /d^xd^. {AW, = (95) 

These extrema are given by solutions of 

2e'^dj{Ba/p) =0 => di{Ba/p) = and d2iBa/p) = (96) 



dA, 



The only solutions are B = (cp/a) = (cg/p) where c is a constant and g = det gij. These 
extrema lie on leaves where the Casimirs take the values 

^^"/i = 0"! a-"" pd^x (97) 

We do not yet understand the physical meaning of these extrema of energy. They have a finite 
energy if i? = {cp/a) vanishes at infinity sufficiently fast. 

^Finding B(r) for given 7„,p(r) is similar to the Classical Moment Problem. 
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8 Discussion 



A summary of the paper was given in the introduction Sec^ Here, we mention a few directions 
for further study. We would hke to know whether there is a deeper integrable structure that 
would explain the presence of an infinite number of conserved charges in involution for the 
uniform volume measure ^ . The extension to an arbitrary volume form seems likely. Are there 
any time-dependent exact solutions of the non-linear evolution equation? What is the Poisson 
algebra of loop observables and can the Hamiltonian be written in terms of them? Can this 
model be quantized and is there a non-abelian extension? An extension to 3 -|- 1 dimensions is 
possible, though the Poisson algebra has only one analytically known Casimir, the Hopf or link 
invariant. Is there a Lorentz covariant theory along these lines? Can the idea that gauge fields 
be thought of as dual to volume preserving vector fields be exploited in any other context? 

How is our gauge theory related to hydrodynamics and turbulence? Recently, Jackiw, et. 
qj.pH] have studied perfect fluids and certain non-abelian extensions. Our gauge theory shares 
the same phase space as ideal hydrodynamics, but the two theories have different hamiltonians. 
However, similar methods may be useful in the study of both theories. For example, Iyer and 
Rajeevj^ (see also Sec. 11. D of Ref. jTJ]) have proposed a statistical approach to 2 dimensional 
turbulence, based on a matrix regularization of the phase space. It may be possible to use a 
similar regularization for our gauge theory. For an N x N matrix regularization to be integrable, 
it would appear that we need 0{N'^) conserved quantities, while /„ and Hm furnish only 2A^ 
conserved quantities. It is unclear what this implies for the integrability of the continuum 
theory we have proposed in this paper. On the other hand, Polyakov has proposed a theory of 
turbulence in 2 + 1 dimensions based on conformal invariancefSO]. Since the group of conformal 
transformations and area preserving transformations are disjoint except for isometries, it appears 
unlikely that there is any direct relation of our work to Polyakov's. 
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A Poisson Manifolds and Coadjoint Orbits 

We collect a few facts about classical mechanics that we use, to make the paper self- 
contained and fix notation. 

A.l Poisson Manifolds and Symplectic Leaves 

The basic playground of classical mechanics is a Poisson manifold. It is a manifold M with 
a product {., .} : !F[M) x J^{M) — > J^{M) (the Poisson bracket p.b.) on the algebra of 
observables. {., .} is bilinear, skew symmetric and satisfies the Jacobi identity and Leibnitz rule. 
J-{M) is a class of real- valued functions (say, C°°{M)). Any such function / : M — > R generates 
canonical transformations on M . A canonical transformation is a flow on M , associated to the 
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canonical vector field Vf . The Lie derivative of any function g{A) along the flow is given by 
the p.b. Cvfg{A) = {f,g}{A), A G M. Flow lines of the canonical transformation generated by 
/ are integral curves of Vf . 

Often, Poisson algebras of observables are degenerate. They have a center (Casimirs) which 
have zero p.b. with all observables. In such a situation, the Poisson manifold as a whole cannot 

serve as the phase space of a physical system, since it would not be a symplectic manifold. 
Rather, it is the symplectic leaves of a Poisson manifold that can serve as phase spaces. The 
symplectic leaf of a point ^ G M is the set of all points of M reachable from A along integral 
curves of canonical vector fields. On a symplectic leaf, the Poisson structure is non-degenerate 
and can be inverted to define a symplectic structure, a non-degenerate closed 2-form uj . Indeed, 
if ^,,1] arc tangent vectors at ^ to a symplectic leaf, then the symplectic form at A is uj{^,ri) = 
{f,g}{A) where / and g are any two functions whose canonical vector fields at A coincide with 
^ and rj; i.e. ^ = Vf\A,r) = Vg\A- Moreover, on any symplectic leaf, a;(V/,.) = df{.) where df 
is the exterior derivative of /. 

The hamiltonian vector field Vh is the canonical vector field of the hamiltonian H : M ^ Ti. 
Vh generates time evolution ^ = {H,f} = Cv^f ■ Irrespective of the hamiltonian, time 
evolution always stays on the same symplectic leaf. 

A. 2 Coadjoint Orbits in the Dual of a Lie Algebra 

A natural example of a Poisson manifold, that occurs in many areas of physics, is the dual of a 
Lie algebra. The symplectic leaves of the dual of a Lie algebra are the coadjoint orbits of the 
group. To understand this, suppose G is a group, Q its Lie algebra and Q* the dual of the 
Lie algebra. Then we have a bilinear pairing between dual spaces (A, n) G R for A £ Q* and 
u € Q . Suppose f,g are two real-valued functions on Q* . Then their p.b. is defined using the 
differential df{A) G G 

{f,g}{A) = {A,[df,dg]) (98) 

[df,dg] is the commutator in Q. This turns G* into a Poisson manifold, which is often degener- 
ate. The symplectic leaves are coadjoint orbits of the action of G on Q* . To see this, first define 
the inner automorphism Ag : G ^ G, Agh = ghg^^ which takes the group identity e to itself. 
The group adjoint representation Adg : Q ^ Q \s the linearization of the inner automorphism 
at e: Adg = AgJ^^, and is AdgU = gug~^ for matrix groups. The group coadjoint representation 
Ad*g :g*^g* is defined as {Ad*gA, u) = {A, Adgu) . The group coadjoint orbit oi A e G* is 

Oa = {Ad*gA \geG} (99) 

The Lie algebra adjoint representation is ad„ : Q ^ Q where adu = ^\t=oAdg(j^^ for a curve g{t) 
on the group with g{Q) = e and g{Q) = u and takes the form ad^v = [n, v] for matrix groups. 
The Lie algebra coadjoint representation ad* : Q* ^ Q* is defined by [ad^^A, v) = (A, aduv) . 
The Lie algebra coadjoint orbit of A G G* is the tangent space at A to the group coadjoint orbit 
of A. 

{aduA I n G a} = TaOa (100) 

Thus, a tangent vector to a coadjoint orbit at A may be written as ^ = ad^^A for some (not 
necessarily unique) u G G . 
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There is a natural symplectic structure on coadjoint orbits, which turns them into homo- 
geneous symplectic leaves of Q* . The symplectic form (Kirillov form) lv acting on a pair of 
tangent vectors to the orbit Oa at A is given by io{^,r]) = {A, [u,v]) where u,v G G are any 
two Lie algebra elements such that ^ = ad^^A and rj = ad^A. Thus, the coadjoint orbits are 
symplectic manifolds which foliate G* in such a way as to recover the Poisson structure on the 
whole of Q* . The different coadjoint orbits in Q* are not necessarily of the same dimension, but 
are always even dimensional if their dimension is finite. 

The canonical vector field Vj of an observable / : ^* ^ R at a point A & Q* is given by 
the Lie algebra coadjoint action of the differential df{A), Vf{A) = ad'^j^A 

CvMA) = {/, g}{A) = {A, [df, dg]) = {A, ad^fdg) = {ad*afA, dg) (101) 

In particular, infinitesimal time evolution is just the coadjoint action of the Lie algebra element 
dH , the diff'erential of the hamiltonian. For example, if fiui^) = i^^u) is the moment map for 
u £ G , then the canonical vector field V^^ (A) = ad^^A . In other words, infinitesimal canonical 
transformations generated by moment maps are the same as Lie algebra coadjoint actions. 

Observables in the center of the Poisson algebra (Casimirs) are constant on coadjoint orbits. 
They are invariant under the group and Lie algebra coadjoint actions. To show that an observable 
is a Casimir, it suffices to check that it commutes with the moment maps which generate the 
coadjoint action. 

Example: Eulerian Rigid Body Let il* be the components of angular velocity of a rigid body 
in the co-rotating frame. The components of angular velocity lie in the Lie algebra Q of the 
rotation group G = 50(3). The dual space to angular velocities consists of angular momenta 
Li with the pairing (or moment map) (L, O) = LjO*. The space of angular momenta is the 
dual S0{3)* = R^. The latter carries a Poisson structure {Li,Lj} = eij^L^. Observables are 
real-valued functions of angular momentum f{L) and satisfy the p.b. 

{/.,}(L)=i:{L.,L,}^^ (102) 

The space of angular momenta R"^ must be a degenerate Poisson manifold, since it is not 
even dimensional. Indeed, the symplectic leaves are concentric spheres centered at the point 
Li = as well as the point Li = 0. These symplectic leaves are the coadjoint orbits of 
5*0(3) acting on the space R^ of angular momenta. The symplectic form on a sphere of radius 
r is given by r sin 9d6 A d(p. The Casimirs are functions of L^ = Y^iLl and are constant 
on the symplectic leaves, they are invariant under the coadjoint action of the rotation group 
on R^. The hamiltonian is H = J^i^f/^h if the axes are chosen along the principle axes 
of inertia. li are the principle moments of inertia, the eigenvalues of the inertia operator 
lij '■ Q ^ Q* which maps angular velocities to angular momenta Lj = lijVL^ . The equations of 
motion L = ad*^jjL = ad*^_i^L = adf^L = L x J7 are Lj = {H,Li}, Li = 023.^2-^3 and cyclic 
permutations thereof, where aij = IJ^ — I^^ . 



B The Charges In are in Involution 



We show that the charges /„ are in involution /„} = 0. To calculate 



JM 



Sim Sin 



5A 



pSA, 



Sin A 1 Sim 



5A 



pSAi 



d^x 



(103) 
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sir, 



we need 

So the p.b. becomes 

{/„,/„} = mne'^e^^ [ (fx Ai 
Jm 



me'^dj{{B/p)"'-') 



1 



dk{B/pr-'dA-di{B/p) 



m <^ n 



(104) 



(105) 



The first term in square brackets dk{B / p)"^ ^dj\^di{B / p)'"' "^j can be written as 

_ i)(„ _ l){B/pr-^dk{B/p)d,\^p-\B/pT-^di{B/p) 

= (m- l){n - l){B/pr-^^-%{p-^)dk{B/p)di{B/p) 

+ (^ _ i)(„ _ i)(n _ 2){B/pr+^-'>p-'dk{B/p)dj{B/p)di{B/p) 

+ (^ _ i)(„ _ l){B/pr+^-^p-'dk{B/p)djdi{B/p) (106) 

The I''* and 3'''^ terms are symmetric mider in ^ n and therefore do not contribute to {Im, In} ■ 
Therefore we get (/i = 0, so we can ignore m,n= 1) 



{Irrtj In} 



mn{m — l)(n — l)(n — m) 



m+n—5 



e''e^''AMB/p)d,{B/p)dk{B/p) 



Now the term in square brackets vanishes identically due to antisymmetry of e^^ . We conclude 
that {/mi In} = 0. Thus, In, n = 1,2,3 .. . are an infinite number of charges in involution. 



C /2 is a Casimir of the Poisson Algebra 

To find the infinitesimal change of In under the coadjoint action, we calculate {In,Pu}- If this 
vanishes for all volume preserving u, then /„ would be constant on symplectic leaves and hence 
a Casimir. Recall H28() that the differential of 7„ is 



n-l 



{dIn{A)y = np-^e'^d,{B/p) 



and the differential of the moment map is {dp^Y = ■ Their Lie bracket is 

[u,dInX = u^dj{np-h"'dk{B/pT~^)-np-h^^dk{B/pr-^d,u\ 
= ny^uidj{p'^dk{B/pY-^) - p-^e^^dk{B/pY-^djU^ 



Thus {pu,In} = J Ai[u,dInYpd^x gives 



{fJ'u,In} = n J d^xA. 
For p = 1 this becomes 

{lJ'u,In} =n d^xAi 



^^pu^dj{p-^dk{B/pY-^) - e^'^dkiB/py-^djU^ 



Y'^u^idkdjB''-^) - £^\dy)idkB"-^) 



(107) 



(108) 



(109) 



(110) 
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Specializing to n = 2 , and writing u in terms of its stream function, = £*'5/'0 , 
{/2,M = 2 J (fx e'^e^''Ai{didj^lj)idkB) - £'^e^'AiidkdjB)idi4>) 



(111) 



{-^2;/Uu} is gauge-invariant, so we calculate it in the gauge Ai = 0, B = diA2 and denote 
A2 = A, = X, x^ = y and derivatives by subscripts. 



^{I2,l^u} = J dxdyA 



'4^yAxxx '4^x-^xxy ~l" '4^xyAxx '4'xx-^xy 



(112) 



The idea is to integrate by parts and show that this expression vanishes. Let us temporarily call 
the second factor of A by the name A. Write |{/2,//u} as a sum of four terms T1 + T2 + T3 + T4 



^xxy 



Ti = j dxdyAt/jyAxxx; T2 = - j dxdyAtpxAx 
^3 = j dxdyAx/jxyAxx; T4 = - j dxdyAtpxxAxy (113) 

We will show that Ti +T2 = —T^ — T4. Integrating by parts till there are no derivatives on A, 

Ti = —J dxdyA^^yAxxx ~^ '^Axx'4^xy '^-^x4'xxy -^"^^xxxy^ 

T2 = J dxdyA^^xAxxy + "^Axytpxx + "^Axtpxxy + Axx'fpxy + Aytpxxx + Atpxxxy^ 

T3 = J dxdyA(^^pxyAxx + ^Axtpxxy + All^xxxy^ 

T4 = —J dxdyA^Axylpxx ^y'i'xxx '''' ^x'i'xxy ^'4^xxxy^ 0-^^) 

Using the fact that A = ^4 we get for Ti and T2 

J dxdyA^Axx'fpxy + ^■^x'^Pxxy + ■^''Pxxxy^ 
2 X r2 = j dxdyA^Axytpxx '^^x'4'xxy -^xx'^xy ■^y'4^xxx -^"^xxxy^ (115) 

Ts and T4 give us the identities 

J dxdy AAxi'xxy = - J dxdy AA^xxxy 

j dxdy AAyijjxxx = j dxdy AAxipxxy (116) 
Use these to simplify Ti and T2 by eliminating ipxxxy smd ipxxx in favor of tpxxy '■ 
2 X Ti = - J dxdy A(^Axxipxy + Axtjjxxy^ 

2.T2 = jdxdyA{2Axy^xx+Axxi^xy + Ax.^xxy) (117) 
Adding these and setting A = A 

n+T. = J^iyA[A..,^,.-A,,i,..) (118) 
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Meanwhile by definition, 

Thus {l2,l^u} = 2{Ti+T2 + n + n) = (119) 
We conclude that I2 lies in the center of the Poisson algebra for constant p . 

D H = ^ f(dA/iiyaii is not a Casimir 

I2 = J{dA/ fj,)^fx turned out to be in the center of the Poisson algebra of gauge-invariant 
functions. Here we show that the magnetic energy H = ^ J {dA / p,)^ a p, with a non-constant 
a = p^/g, does not lie in the center, and therefore leads to non-trivial time evolution. We do 
this by giving an explicit example of a gauge-invariant function with which it has a non- vanishing 
p.b. Consider p = I, then from (|55|) 



{H,pu} = j d'^xAi e^^{djU^)dk{Ba)-e'^u^{djdkBa) 
In gauge Ai = 0, A2 ^ A and with = e'^Wjip we get {x^ = x,x'^ = y 



(120) 



{H,pu} = J dxdyA 



+{dyi^)dl{Ba) - id,^lj)d,y{Ba) 



:i21) 



Now for the simple choices 'ip = xy,A = a = e +y we have B = dxA = —xe 
The p.b. can be calculated exactly to yield 

{H, p^} = J dxdye-'^^'+y'y' ^-l-Ax'-2y' + Ax\2 + y^)) = ^ (122) 

Thus H does not lie in the center of the Poisson algebra. We also checked that H transforms 
non-trivially under many other generators pu of the coadjoint action. 



E Hn are in Involution for Uniform Measure 

Suppose p is the uniform measure {p = 1). We prove that the charges 

Hn = j{dA/pYap = j {B/pY(Tpd^x = J B"ad'^x (123) 

are in involution 

{H^,Hn} = for m,n = 0,l,2,3,... (124) 

An immediate corollary is that Hn are conserved quantities since the hamiltonian is ^H2. Here 
S ^ at 00 and a = (p/Qg)'^ = 1/5 must be such that these integrals converge. The proof 
involves explicitly computing the p.b. and integrating by parts several times. Using 

^ = 5^Hn = ne'^d,{a{B/pr-') (125) 
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we can express the p.b. as 



X 

For p = 1 this becomes 

(m + l)(n + 1) 



A. 



(fx = mne^^e^^ / A 



dk{a{B/pr-^)di{p-^di{a{B/pr-^)) -m^n 



(fx 



A, 



dk{aB"')djdi{aB'') -m^n 



(fx 



(126) 



(127) 



Now we expand out the derivatives of products of B and a and ehminate terms that vanish 
due to antisymmetry of e^^ or antisymmetry in m and n . We get 



mn{n — m) 



B{dka){dia){djB) + aB{dka){djdiB) 



+{n + m- 3)a{dka){djB){diB) + aB{djB){dkdia) 



(128) 



Since this p.b. is gauge-invariant, calculate in the gauge Ai = 0, call A2 = A, B = d^A and 
denote derivatives by subscripts {x^ = x,x'^ = y) 



{Hm+i,Hn+i} = -e^^m + l)(n + l)(n - m) j dxdy A5'"+"-2 Ba^B^ak 
+uBBjxUk + (n + m — l)aBxBjak + uBBjUkx 
After collecting terms, this becomes 



(129) 



{H^, Hn} = mn{m-n) J ^5"*+""^ 



{B^ay - Bya,^){Bax + (n + m - 2>)Bxa) 



"^CfB-I Bxx^y -^xy^x ~l~ -^x^xy -^y^xx 



(130) 



To simplify it we define k = m + n — A and use identities such as B'^B^ = {B^^^)^/ {k + 1) and 
(^Cx = |(c^)x to combine the factors of B and a to write 



{Hm,Hn} = 



mn{m — n) 

T3 = i I A{B^+%B,{a\ 
n = \j AB^+^B,,{a\ 



Tj where 

i=l 
T2 



T7 



k + : 



xO'CTxy 



n = -\j A{B'^+\Bx{a% 
^ I AB^+^Bxy{a% 

\-2 j ^(^'^') 



Is 



k + ' 



y(T(Txx 



(131) 



Numerically we find Ti + + T2 + Tg = —{T^ + Ts + T4 + Tg) . To see this analytically, integrate 
by parts with the aim of eliminating A in favor of 5 = dxA . Using similar identities as before, 



1 /• (S'=+3),(a^), 1 f A{B>'+%x(.a') 



k + 3 



k+2\ 



k + 2 
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T2 = -n + 



1 ( iB'^^%ia% 1 f 

2 J k + 3 2 J k + 2 

^ ^ 2j A; + 3 2j k + 2 

^ 2j k+2 2j k+2 ^ ' 



Then 



T^ + n + TA+ Tq 



1 f {B>^+%ia')y 1 f AyiB'^+%ia% 
2J A; + 3 2J k + 2 

T^ + Tj + T^ + T, = j Ai^{B'^^Xy{a% - {B^+%.{a\ 

{Hm,Hn} = mn{m-n){Ui + U2 + U3 + U4) (133) 

where + + T4 + Tq = Ui + U2 and Ti + Tt + T2 + Tg = U3 + U4. Finally, integration by 
parts shows that 

and that U2+U3 = —U4—U1. We conclude that J^f^i Ui = and therefore {Hm,Hn} = 0. Since 
the hamiltonian is half the second charge, H = \H2 , we have shown that Hn for n = 1, 2, 3, . . . 
are an infinite number of conserved quantities, which moreover, are in involution! 
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